Interference of an array of independent Bose-Einstein condensates 
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We have observed high-contrast matter wave interference between 30 Bose-Einstein condensates 
with uncorrelated phases. Interference patterns were observed after independent condensates were 
released from a one-dimensional optical lattice and allowed to expand and overlap. This initially 
surprising phenomenon is explained with a simple theoretical model which generalizes the analysis 
of the interference of two independent condensates. 
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PACS numbers: 03.75.Lm, 32.80.Pj 

Studies of Bose-Einstein condensates (BECs) loaded 
into the periodic potential of an optical lattice have been 
continuously growing in the recent years 0. These sys- 
tems have a great potential for a range of applications 
such as modelling of solid state systems 0, IE @ , 
preparation of low-dimensional quantum gases [7|,|8|, tun- 
ing of atom properties 0, El j trapped atom interferom- 
etry and quantum information processing 

One of the most commonly used probes of these novel 
systems are the interference patterns obtained when the 
gas is released from the lattice, so that the wave packets 
emanating from different lattice sites expand and over- 
lap In particular, the appearance of high-contrast 
interference fringes in the resulting density distribution 
is commonly associated with the presence of phase co- 
herence between different lattice sites. 

In this Letter, we study the interference of a regular 
array of Bose-Einstein condensates with random relative 
phases. Independent condensates, each containing ~ 10 4 
atoms, are produced in a one-dimensional (ID) optical 
lattice, in the regime where the rate of tunnelling between 
the lattice sites is negligible on the time scale of the ex- 
periment. Contrary to what could be naively expected, 
we show that high-contrast interference fringes are com- 
monly observed in this system. We present a theoretical 
model which quantitatively reproduces our experimental 
results, and show that the periodicity of the lattice is 
sufficient for the emergence of high-contrast interference 
patterns, even in the absence of phase coherence between 
the condensates. This conclusion is independent of the 
number of sites or the dimensionality of the lattice. Our 
results generalize the analysis of the interference between 
two independent condensates (lilll6lll7| . In the last part 
of the paper we briefly discuss the potential of our system 
for creating strongly number-squeezed states and report 
on an unexplained heating effect which occurs for a nar- 
row range of lattice depths. 

Our experiments start with a quasipure 87 Rb conden- 
sate with 3 x 10 5 atoms in the F = nip = 2 hyper- 
fine state. Condensates are produced by radio-frequency 
evaporation in a cylindrically symmetric magnetic trap. 
The harmonic trapping frequencies are lo±/2tt = 74 Hz 
radially, and ui * 1 /2n = 11 Hz axially, leading to cigar- 




FIG. 1: Interference of Bose-Einstein condensates with un- 
correlated phases, (a) A deep ID optical lattice splits a cigar- 
shaped condensate into 30 independent BECs. (b) Absorption 
image of the cloud after 22 ms of expansion from the lattice. 
The density distribution shows clear interference fringes of 
high contrast, (c) Axial density profile of the cloud, radially 
averaged over the central 25 /im. The fit described in the text 
gives a fringe amplitude A\ = 0.64 for this image. 

shaped condensates with a Thomas-Fermi length of 
Ltf = 84 ^m, and a radius of i?TF = 6 /im. 

We then ramp up the periodic potential created by a 
ID optical lattice. The lattice is superimposed on the 
magnetic trapping potential along the long axis (z) of 
the cigar (FigGJi). Two equally polarized laser beams of 
wavelength A = 532 nm intersect at an angle 8 = 0.20 rad 
to create a standing wave optical dipole potential with a 
period of d — A/2 sin(0/2) = 2.7/im. The two beams 
are focused to waists of about 100 /im, and carry a laser 
power of up to 220 mW each. The blue-detuned laser 
light creates a repulsive potential for the atoms which 
accumulate at the nodes of the standing wave, with the 
radial confinement being provided by the magnetic po- 
tential. Along z, the lattice potential has the shape 



Vn 

V{z) = y cos{2n z/d) 



(1) 



The lattice depth at full laser power is Vo/h ss 50 kHz, 
and the number of occupied sites is N = Ltf /d ~ 30. 

Thanks to the long lattice period and the large hight 
of the potential barrier between the sites, the 30 conden- 
sates can be completely isolated from each other. The 
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FIG. 2: Polar plots of the fringe amplitudes and phases 
(Ai,Bi) for 200 images obtained under same experimental 
conditions, (a) Phase-uncorrelated condensates, (b) Phase- 
correlated condensates. Insets: Axial density profiles aver- 
aged over the 200 images. 

matrix element J for the tunnelling between neighboring 
sites scales as E R e^ 2 ^ Vo/ER , where E R = h 2 k 2 /2m = 
h x 80 Hz, k = Tr/d, and m is the atom mass. Our max- 
imum lattice depth corresponds to Vq ~ 600 Er, leading 
to a timescale h/ J of tens of years. 

At full lattice depth, the gas at each site is in a quasi 
2D regime. The motion along z is "frozen out" because 
the oscillation frequency ui z /(2ir) is 4kHz, while the tem- 
perature and the chemical potential of the gas correspond 
to frequencies smaller than 2.5 kHz. Each condensate is 
therefore in the harmonic oscillator ground state along 
this direction, with a density distribution given by a 
Gaussian of width I — y/H/ (2mw z ) = 120 nm <C d. 

In a typical experiment, we ramp up the lattice to 
600 Er in r r amp = 200 ms. After holding the atoms in 
the lattice for Thoid = 500 ms, the optical and magnetic 
trapping potential are switched off simultaneously, and 
the density distribution of the cloud is recorded by ab- 
sorption imaging after t = 22 ms of timc-of-flight (ToF) 
expansion. The hold time Thoid is sufficient for the phases 
of the independently evolving condensates to completely 
decorrelate (see below). Despite this, the images com- 
monly show clear interference fringes. A spectacular ex- 
ample of this surprising phenomenon is given in Fig.^>-c, 
where the contrast is > 60 %. Note that high-contrast 
interference is also observed if the lattice is ramped up 
before the end of the evaporation sequence, so that the 
condensates are produced independently at different lat- 
tice sites and "have never seen one another" 

We analyze the images by fitting the axial density pro- 
files (Fig. [Tt) with [1 + Ai cos(Si + 2-kz/D)} G{z), where 
G(z) is a Gaussian envelope. This procedure extracts 
the first-harmonic modulation of the density distribution, 
with the fitted fringe period of 39 /im in agreement with 
the expected value D = ht/(md). In Fig. [3]}, we sum- 
marize our results for 200 consecutively taken images. In 
most cases the fringe amplitude A\ is significant, with 
the mean value (Ai) — 0.34, and the standard deviation 
a A t = 0.17. The fringe phase B\ is randomly distributed 
between and 27r. Consequently, no periodic modulation 



remains visible if we average the 200 density profiles. 

In Fig.Eb, we contrast this with the interference of 30 
phase-correlated condensates. Here, we ramp up the lat- 
tice in Tramp = 3 ms and immediately release the BECs 
from the trap, before their phases completely diffuse away 
from each other. In this case, the fringe phases B\ are 
clearly not randomly distributed, and the sum of 200 im- 
ages still shows interference fringes of pronounced con- 
trast. Phase correlations between the separated conden- 
sates were lost if the lattice was left on for Thoid _ 3 ms. 

We have numerically simulated our experiments 
with phase-uncorrelated condensates using the following 
model. We consider a ID array of N = 30 BECs ini- 
tially localized at positions z n = nd (n = 1, . . . , N), and 
assume that each condensate is in a coherent state de- 
scribed by an amplitude a n and a phase <j> n 0] . We set 
a n oc n(N — n), corresponding to the Thomas-Fermi pro- 
file of our BEC at the point when the lattice is switched 
on. For expansion times t ^> l/w z ~ 40 /is, the atom 
density at position z is given by [8|: 



I{z) (X 



N 

E 

n=l 



a n e 



e im(z-z n ) 2 /{2Ht) & - 



(z-z n f/Z% 



(2) 



where Zq = ht/(m£). We neglected the effects of the 
atomic interactions during the expansion. We perform a 
Monte-Carlo analysis of I(z) by assigning sets of random 
numbers to the phases {<j> n }- We convolve the resulting 
density profiles with a Gaussian of 5 /im width to account 
for the finite resolution of our imaging system, and then 
fit them in the same way as the experimental data. The 
fitted fringe phase B\ is randomly distributed between 
and 27T. For the fringe amplitude we find (Ai) slm = 
0.31 and a s F^ = 0.16, in excellent agreement with the 
experiment |20j. 

In order to give a more intuitive explanation of our 
observations, and generalize our analysis to arbitrarily 
large N, we now make the following simplifications. We 
assume that all condensates contain the same average 
number of atoms, a n — a, and that the expansion time t 
is large, so that £, z n <C Z , yjht/m. In this case, Eq. @ 
can be rewritten in the form I(z) oc N a 2 e~ 2z * ' z o F(z), 
where 



JV-l 



F(z) = 1+^2 A n cos(B n + 2-nnz/D) 



(3) 



is a periodic function with period D = ht/(md) and av- 
erage value 1. F(z) contains all the information about 
the contrast of the interference pattern. The ampli- 
tude A n and the phase B n of the n-th harmonic of 
F(z) are given by the modulus and the argument of 

(2/iV)Ef=„ + ie^'-^-»). 

If the N condensates have the same phase, F(z) cor- 
responds to the usual function describing the diffraction 
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of a coherent wave on a grating 



F(z) 



1 sin 2 (Nttz/D) 
N sin 2 (tt z/D) 



(4) 



In this case, F(z) has sharp peaks of hight N and width 
~ 1/N at positions z = pD, where p is an integer. In 
between the peaks, F(z) ~ 0. 

For the case of uncorrelated phases {</>«}, using Monte- 
Carlo sampling, we find for 10 < N < 10 4 : 



(Fmax) ~ 1.2 ln(AT) » 1 (F min ) 



0.2 



< i 



(5) 



where -F m ax and F m i n are the maximum and the mini- 
mum "single-shot" values of F(z), i.e. for a given set 
{<j>n}- This shows that even for very large N, most single- 
shot F(z) are strongly modulated, with a contrast of al- 
most 100%. Since F(z) is also periodic with period D, 
each single-shot has the qualitative appearance of a high- 
contrast interference pattern normally expected in the 
coherent case. However, the exact shape of F(z) and the 
position of its extrema vary randomly from shot to shot. 

The harmonic content of single-shot F(z) is given by 
C n = {A 2 n ) = 4(N — n)/N 2 , each A n resulting from sum- 
ming N — n complex numbers with random phases and 
moduli 2/N. The large modulation of F(z) for N — > oo 
can be understood by noting that while the weight of each 
harmonic decreases, their number increases |2l|. It is in- 
teresting to contrast this with averaging N shots of two- 
condensate interference with an irreproducible phase ■ 
In the latter case, the weight of the first harmonic is simi- 
lar (~ l/ViV), but the absence of higher harmonics leads 
to a vanishing contrast for large N. 

We obtain similar results if we choose for the initial 
state a Mott insulator with exactly one atom per lattice 
site (see also [23J). In this case, the signal after ToF 
consists of the set of coordinates £i , . . . , where the N 
atoms are detected, and its harmonic content is given by: 



JV 



C n 



JV(JV-l)^ 



,i2ir(Ci-C,/W-D 



where (...) now denotes a quantum average. A simple 
calculation gives C n = 4(JV - n)/[N(N - 1)]. 

All our arguments naturally extend into two and three 
dimensions. In 3D, the function generalizing Eq. J2J is 
F(r) = 1 + '}2 n A T iCos{B rl + 2irn ■ r/D) where n — 
(n x ,n y ,n z ) is a triplet of integers. This is a periodic 
function in x, y, z with a period D. A Monte-Carlo anal- 
ysis for a cubic lattice with 20 x 20 x 20 sites with random 
phases shows that F(r) is again strongly modulated, with 
(F max ) ~ 12 and (F mia ) < 1. 

In the experimental study of the superfluid to Mott in- 
sulator transition with cold atoms |3J, |]| i the disappear- 
ance of interference fringes is observed in the insulating 



domain, and used as one of the signatures for the loss 
of long range phase coherence. This seems in contradic- 
tion with our results. However one can reconcile the two 
findings by noticing that the experiments 0, 0] are per- 
formed with a 3D system, observed by integrating the 
spatial distribution I(r) along the line of observation z. 
In this case, only the harmonics (n x , n y , 0) are observed. 
In a 3D experiment performed with N x x JVj, x N z sites, 
the integration along z thus reduces the modulation am- 
plitude by \/N z . In fact, the effect is similar to that of 
averaging N z i mag es in a 2D experiment, performed with 
N x x N y sites 24]. 

In the next part of this Letter, we briefly discuss the 
potential of our setup for producing "number-squeezed" 
states with large occupation numbers 0], where the 
atom number on each site has a sub-Poissonian distri- 
bution with an average value no > 1 and a standard 
deviation a < y / no~. Such states are important for atom 
interferometry and precision measurements |27| . 

For simplicity, we restrict our discussion to the case 
of translational invariance, without the quadratic mag- 
netic potential applied along the z-axis. In the ground 
state of the system, the squeezing of the atom number 
on each site depends on the ratio of the tunnelling rate 
t = noJ/h and the effective strength of the repulsive 
on-site interactions U ~ m/^Oj where /i is the chemical 
potential. To give a sense of scale, in our case no ~ 10 4 
and U/h ~ 0.2 Hz. We can qualitatively distinguish three 
regimes 0IH: (i) For J > no U, squeezing is negligi- 
ble and the atom number on each site follows a Pois- 
sonian distribution (a — ^Juq ~ 100). (ii) For J ~ U 
(t ~ 200 Hz), squeezing is significant and a is reduced to 
n^ 4 ~ 10. (iii) Finally, a phase transition to a Mott insu- 
lator state with a < 1 occurs for J ~ U /tiq (t ~ 0.2 Hz). 

In our setup, we can tune the value of J across this 
full range, and the criterion for the Mott transition with 
no ~ 10 4 atoms per site is satisfied for Vq ~ 100 — 150 En. 
However, we point out two practical difficulties which 
arise for such large values of no. First, for the system to 
be in its ground state, all decoherence processes, such as 
particle loss, must have rates lower than t. Second, the 
assumed translational invariance must be insured to a 
sufficient level, so that the potential energies at different 
lattice sites match to better than n^J = ht. While it 
seems difficult to fulfill these criteria for the values of t 
low enough for the Mott transition to occur, the regime 



of strong squeezing, a 



< 



1/4 



should be accessible. 



In exploring the range of lattice depths 0—150 En, we 
have also observed an unexplained heating effect. In the 
region 20 < Vq/Er < 45 we observe strong heating of 
the system which peaks for Vq/Eh ps 25 and 35 (Fig. - 
The peak heating rate (measured from the radial size of 
the cloud after ToF) is ~ 200nK/s as long as the gas is 
at least partially condensed. However, once the conden- 
sate disappears, heating of the thermal cloud becomes 
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FIG. 3: Axial size of the cloud after 22 ms of expansion as a 
function of the lattice depth Vb- Inset: Zoom-in on the region 
where an unexplained heating of the cloud occurs. For all 



data points the lattice was raised in r ri 



200 ms and left 



on for Ti 10 id = 500 ms before ToF. Error bars are statistical. 
Outside the heating region, the slow increase of the axial size 
matches the expected dependence Zo tx l/£ oc Vq^ 4 . 

negligible (< lOnK/s). Outside 20 < V /E R < 45, heat- 
ing due to the lattice is always negligible, independent 
of the lattice depth or the condensed fraction. We could 
not attribute this heating to any trivial technical effect. 
Also, the dynamical instability which occurs at a finite 
relative velocity between the condensate and the lattice 
(|29j and refs. therein), should not be relevant here. 

In conclusion, we have studied a ID periodic array of 
independent condensates prepared in a deep optical lat- 
tice. We have shown that most "single-shot" realizations 
of this system show high-contrast interference patterns 
in a time-of-fiight expansion. We have explained this ef- 
fect with a simple model which naturally extends to 3D 
lattices. This initially surprising result should be taken 
into account in the ongoing studies of atomic superfluid- 
ity and coherence in optical lattices, where the contrast 
of the interference patterns is often used as a diagnostic. 
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